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I. INTRODUCTION 

We report the results of a search for high-mass reso- 
nances decaying to ZZ in pp collisions at y/s = 1.96 TeV 
at the Tevatron. Although the decay of the standard 
model Higgs boson to ZZ is expected to be beyond the 
sensitivity of the Tevatron experiments Q] , new physics 
could affect ZZ production in different ways. In mod- 
els containing large extra dimensions the ZZ production 
cross section is increased through loop corrections [2J. 
Resonances appearing at high mass such as a Randall- 
Sundrum (RS) graviton [3] could decay manifestly to 
two Z bosons. The original RS model predicts Kaluza- 
Klein excitations of the graviton (G*) that decay predom- 
inantly to a pair of charged leptons or a pair of photons. 
Experimental searches for such high-mass resonance de- 
cays have excluded RS graviton states up to a mass of 
around 1 TeV/ c 2 at 95% confidence level for a natural 
choice of coupling parameter [I], both at the Tevatron 
and at the LHC [5]. However, in RS models that have 
standard model fields propagating in the bulk, the G* 
couplings to light fermions and photons may be heavily 
suppressed so that the dominant decay modes are to tt, 
Higgs pairs, or pairs of heavy bosons [6]. Furthermore, 
in some models the decay to heavy bosons is dominant 
[TJ. Suppression of the couplings to light fermions also 
results in gluon fusion becoming the primary production 
process. 

The CDF experiment has previously searched for res- 
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onances decaying to Z pairs and excluded RS gravi- 
tons with mass up to around 0.5TeV/c 2 at 95% con- 
fidence level [8]. The search described in this paper 
gives improved sensitivity over the previous analysis 
through modified event selection, the inclusion of extra 
final states, and the addition of more data. Three fi- 
nal states are examined, corresponding to the different 
Z boson decay modes ZZ -> £+£-£+£-, ZZ -> £ + £~uv, 
and ZZ — > £ + £~~ jj, where £ is an electron or muon and 
j is a hadronic jet. These three channels have different 
signal-to-background ratios and allow an overconstrained 
search. The four-lepton final state has the smallest back- 
ground; however, depending on the resonance mass, the 
best single-channel sensitivity is provided by either the 
ZZ -> £ + £~jj or ZZ -> £+rvu channels. 

The paper is organized as follows: in Section [IT] we 
introduce the CDF detector and trigger system; in Sec- 
tion |TlT] we describe the reconstruction and identification 



procedures; then in Sections IV VI we report the search 
results from each of the channels ZZ — > £ + £~£ + £~ , 
ZZ -> £+£-vv and ZZ -> £ + £~jj. Section[vTT|gives lim- 



ity is measured by a hodoscopic system of Chercnkov 
counters [T5] . 

CDF has a three-level online trigger system. The data 
used in this measurement were collected using inclusive 
high-p^ electron and muon triggers, and a two-electron 
trigger. The single-lepton triggers select events that have 
electron or muon candidates with px > 18GeV/c and 
\rj\ < 1.0 16 , and the two-electron trigger uses only 
calorimeter information and allows electron candidates 
above the same px threshold anywhere in the detector. 
The data correspond to an integrated luminosity of 6 fb^ 1 
collected between February 2002 and February 2010. 



III. RECONSTRUCTION AND 
IDENTIFICATION 

In this section we discuss lepton reconstruction and 
identification, and reconstruction of jets and missing 
transverse energy. 



its resulting from all three channels and their combina- 
tion. 



A. Leptons 



II. DETECTOR 

The CDF II detector is a general-purpose particle de- 
tector, described in detail elsewhere [5]. The results re- 
ported in this paper use information from several detector 
subsystems for charged lepton and jet reconstruction and 
identification. 

Tracks of charged particles are reconstructed in the sil- 
icon system [TU] and in the central tracker [TT], which is 
a drift chamber that consists of 96 layers of sense wires 
grouped into eight 'superlayers'. Superlayers alternate 
between an axial configuration, with sense wires parallel 
to the colliding beams, and a small-angle stereo config- 
uration. For high momentum tracks the resolution is 
a PT /p T ~ 1.7 x 10~ 3 ( GeV/c)^ 1 , where pr = psin6>, p 
being the track momentum and 9 the polar angle with 
respect to the proton beam direction. 

The calorimeter is segmented radially into electromag- 
netic and hadronic compartments [TSJ, [T3] . The central 
calorimeter is split at the center into two separate bar- 
rels and covers the pseudorapidity range \rj\ < 1.1 (where 
T) = — In tan | ) . Each barrel consists of 24 azimuthal 
wedges segmented in projective towers of 0.1 in 77. The 
forward calorimeter segmentation increases from 0.1 in 
77 and 7.5° in the azimuthal angle <f> at 77 = 1.1, to 0.5 
in r\ and 15° in <f> at r\ = 3.6. Electron energy reso- 
lutions are 13.5%/-\/Et © 2% in the central calorimeter 
and 16%/y/E © 1% in the forward calorimeters, where 
Et = Esm9. The electromagnetic calorimeters incorpo- 
rate shower maximum detectors that are used to measure 
shower profiles with spatial resolution of around 2 mm. 

Dedicated muon detectors [14] are mounted around the 
calorimeters, providing coverage for \rj\ < 1.5. Luminos- 



Decays of a heavy resonance to ZZ, where at least 
one of the Z bosons decays leptonically, result in a wide 
lepton energy spectrum. Any inefficiency in lepton re- 
construction and identification is raised to the fourth 
power in the ZZ —> £ + £~£ + £~ channel. Thus, keeping 
efficiency high while maintaining stringent background 
rejection is equally important for px ~ 20GeV/c and 
for pt > lOOGeV/c. To this end, this analysis incorpo- 
rates several refinements in the offline reconstruction and 
identification of electron and muon candidates. Studies 
were performed on inclusive Z — > £ + £~ candidates and on 
events containing one lepton plus two additional tracks 
having p T > 10 GeV/ c, and this latter data set was fully 
reprocessed for the ZZ — > £ + £~£ + £~ analysis. 

First we describe the elements of the lepton selection 
that are standard to CDF. Electron candidates consist 
of a calorimeter cluster matched to a well-reconstructed 
track. Candidates are required to be within the fidu- 
cial region of the shower maximum detectors and have 
a shower that is mostly contained in the electromag- 
netic compartment of the calorimeter, with a shower 
shape that is consistent with test beam expectation 
[17] . For candidates reconstructed in the central part 
of the detector (|/7| < 1.1), the matched track must 
have pT>WGeV/c, pass through all layers of the cen- 
tral tracker, and have a fit x 2 /d.o.f. < 3. Candi- 
dates reconstructed in the forward part of the detector, 
1.13 < 1 77 1 < 2.8, must either have a track in the central 
tracker, or a track in the silicon system with > 5 hits. 

A muon candidate is reconstructed from a track in the 
central tracker pointing to track segments in the muon 
chambers. Muon track trajectories must be such that 
at least 30 central tracker hits would be expected ge- 
ometrically, and at least 60% of those must be found. 
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Tracks pointing forward that have fewer than three cen- 
tral tracker segments must also have at least five r — <f> 
hits in the silicon tracking system. Muon energy deposi- 
tion must be consistent with that of a minimally-ionizing 
particle. We also consider minimally-ionizing tracks that 
have no track segments in the muon systems as muon 
candidates. 

Electron and muon candidates are required to have 
Et > 15GeV and px > 15GeV/c respectively. In addi- 
tion, one of the lepton candidates in each event is required 
to have E T > 20GeV (electrons) or p T > 20GeV/c 
(muons), and to pass more restrictive quality require- 
ments. These extra requirements are that the lepton 
track must have at least three segments reconstructed in 
the axial superlayers and three in the stereo superlayers; 
and the track of a muon candidate must also be well- 
matched to a track segment reconstructed in the muon 
system. 

The first refinement in lepton selection is in the iso- 
lation requirement made on all lepton candidates. The 
'isolation energy' is the amount of energy reconstructed 
in a cone of A R < 0.4 ar ound a lepton candidate, where 
AR = J (Arj) 2 + (A<j)) 2 . In computing the isolation en- 
ergy, we refine the treatment of energy leakage across 
calorimeter cell boundaries. In the central calorimeter, 
electron clusters include energy depositions from only a 
single wedge in <j>. As each calorimeter tower is read 
out from different <p sides by two photomultiplier tubes, 
the relative heights of the pulses locate the energy de- 
position in <f>. Locating the center of the energy depo- 
sitions in towers neighboring the electron cluster allows 
us to estimate the leakage, and correct the isolation en- 
ergy variable event-by-event, rather than by applying an 
average correction. The correction method is validated 
by examining the isolation energy as a function of shower 
position in the calorimeter cell, which is found to be more 
uniform than under application of the standard average 
correction, as shown in Fig. [lja). Muons are not ex- 
pected to result in energy leakage; their isolation energy 
is also shown in Fig. [lja) as validation of the method. 
The average isolation energy should depend on the in- 
stantaneous luminosity but not on the lepton Et, and 
its uniformity in lepton Et is confirmed by Fig. Qb). 
All electron and muon candidates are therefore required 
to be isolated in the calorimeter by limiting the isolation 
energy to be below 4 GeV. Cutting on isolation energy, 
rather than requiring the standard ratio of isolation en- 
ergy to lepton momentum to be < 0.1 [T7], increases the 
acceptance for ZZ -> £ + £~£ + £~ events by 4%. 

For the ZZ — > £ + £~£ + £~ analysis, events have been re- 
constructed with an updated version of the CDF tracking 
code that gives improved pattern recognition at high lu- 
minosities. The updated version includes an extra algo- 
rithm to associate hits in the central tracker with silicon- 
only tracks from electron candidates in the forward re- 
gion of the detector. Adding extra hits on to these tracks 
improves the robustness of forward electron charge iden- 
tification. 
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FIG. 1: (a) Corrected isolation energy across the calorime- 
ter wedge coordinate X in Z — > e + e~ (new correction: solid 
circles; standard correction: open squares) and Z — > fi + fi~ 
(new correction: open stars) events; (b) average calorimeter 
isolation energy as a function of electron Et in Z — > e + e~ 
events. 



Use of an improved reconstruction algorithm in the 
central shower maximum detector gives better separa- 
tion between showers generated by electron tracks and 
showers produced by bremsstrahlung photons. Matching 
tracks to the showers they initiate in both coordinate and 
energy improves hadron rejection and allows the inclu- 
sion of electron candidates that lose a significant amount 
of energy through bremsstrahlung. The improved back- 
ground rejection allows the relaxation of other standard 
electron identification requirements and, overall, the se- 
lection efficiency is increased by around 9% per electron. 

Electrons reconstructed in the edge 0-rings of the 
calorimeter on cither side of the gap between the cen- 
tral and forward detectors are generally excluded from 
analysis. They are included here, after verification that 
they have energy resolution comparable with electrons 
reconstructed in the bulk of the detectors, and are well- 
modeled in the simulation. This increases electron ac- 
ceptance by around 10% per electron. 

The combined effect of the refinements described above 
is to increase lepton acceptance without increasing fake 
lepton backgrounds, as measured by jet-to-lepton fake 
rates in inclusive jet datasets. The lepton selection 
used for this analysis is validated by measuring inclu- 
sive Z —> £ + £~ cross-sections and separating events by 
calorimeter region and muon system. We verify that for 
each subset of events the measurement is stable in time, 
and combining all channels we measure <j(pp — > Z) x 
Br{Z -> £+£-) = 247 ± 6 (stat.+syst.) ± 15 (lumi.) pb, 
consistent with CDF's measurement 1161. 



B. Jets and fi T 

Jets are reconstructed as clustered energy depositions 
in the calorimeter using a fixed cone clustering algorithm 
with cone size AR = 0.4 [18]. Jet energies are corrected 
for 77-dcpendent calorimeter response and for multiple in- 
teractions P2]. We consider jets having Et > 20 GeV. 

The missing transverse energy (E T ) is defined as the 
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sum over calorimeter tower energies fi T = — X^-^T n i> 
where rii is the unit vector in the transverse plane that 
points to calorimeter tower i. The fi T is adjusted to ac- 
count for the energy corrections made to reconstructed 
jets, and for muons identified in the event. As neutri- 
nos pass through the detector without depositing energy, 
large fi T in an event can imply the presence of high- 
energy neutrinos. 



result of the M«>76 GeV/c 2 requirement, the contribu- 
tion of Z — y tt decays is negligible. 




Z/l* 1 




IV. ZZ -> l+t't+l- CHANNEL 



FIG. 2: Lowest-order standard model ZZ production. 



The first search channel is ZZ £+£-£+£-. We select 
events with four candidate charged leptons, which may be 
electrons or muons. At least two of the four must have 
E T > 20GeV for electron candidates (p T > 20GeV/c 
for muon candidates) and pass the more restrictive lep- 
ton selection; and in order to have the trigger efficiency 
well-defined, at least one must satisfy the trigger require- 
ments. 

Leptons of the same flavor are paired to form Z candi- 
dates, seeded by a lepton that passes the tighter selection. 
In the case of four-electron or four-muon candidates, the 
pairings that minimize the y 2 of the ZZ hypothesis are 
chosen: 

X 2 = (M 12 ~M z ) 2 /a 2 M + (M 34 ~M z ) 2 /a 2 M , 

where M12 and M34 are the masses of the lepton pairs, 
&m = 3 GeV/c 2 approximates experimental resolution 
in Mm for both electron and muon decays, and Mz is 
the mass of the Z boson. 

We find ten events that pass the four-lepton selec- 
tion. In all of these events the number of leptons of 
the same flavor is even. The best pairings of the ten 
candidate events are all oppositely-charged. To mini- 
mize the effect of Z/7* interference, both Z boson can- 
didates are required to be within 15 GeV/c 2 of the Z 
pole, 76 < Mu < 106 GeV/c 2 . Following this require- 
ment, eight event candidates remain: two events have 
four reconstructed electrons (eeee), three have two elec- 
trons and two muons (ee/i/x), and the remaining three 
have four reconstructed muons (fifififi). The two events 
that fail the Z mass requirement both have one Z candi- 
date with invariant mass below 60GeV/c 2 . 

We use the selected events to measure the pp — > ZZ 
production cross section. On- and off-shell ZZ produc- 
tion, as shown in Fig. [5] followed by Z boson decays to 
charged leptons, is the only lowest-order standard model 
process that results in a final state with four high-p^ 
leptons produced in the primary interaction. The back- 
ground in this channel thus comes only from misidenti- 
fication. The main contributions are: pp — > WZ + jet 
with a jet misidentified as a lepton; pp — > Z + 2 jets with 
both jets misidentified as leptons; and pp —> Z + 7 + jet 
with both the photon and the jet misidentified as elec- 
trons. The contribution from tt production is an order 
of magnitude smaller than that of WZ production. As a 



The pythia event generator [20] and the full CDF 
detector simulation |21j are used to simulate kinemat- 
ics of these processes and photon-to-lepton misidentifica- 
tion. Jet-to-lepton misidentification rates are measured 
in inclusive jet data and found to be of the order of 
10~ 4 - 10~ 3 per jet for 15 < E T < 100 GeV. These 
misidentification rates are used to weight the simulated 
events of the background processes, resulting in a total 
background yield estimated to be less than 0.01 event. 

The acceptance for standard model pp — > 
Z/j* Z/Y -» £+£-£+£- is determined using the 
leading-order pythia generator and found to be 
0.17 ± 0.02. The uncertainty has contributions from 
higher-order generator effects, lepton identification, and 
trigger efficiency uncertainty. In order to estimate the 
uncertainty arising from higher-order generator effects, 
the MC@NLO generator [22] is used, interfaced to HERWIG 
[23] to provide parton showering and hadronization. The 
corresponding relative uncertainty on the acceptance is 
estimated to be 2.7%. Lepton identification efficiencies 
are measured in the data using candidate Z — > £ + £~ 
events with uncertainties at the level of 1%. We also ac- 
count for a small drop in lepton identification efficiency 
with time and assign a 2% relative uncertainty per 
lepton for residual run-dependent effects. We assume no 
correlation between the uncertainties on electron and 
muon reconstruction, and full correlation between the 
uncertainties for leptons of the same flavor. The trigger 
efficiency per four-lepton event is close to unity, with a 
systematic uncertainty of less than 0.5%. 

Given the branching fraction for Z — > £ + £~ = (3.366 ± 
0.002)% [24], the branching fraction for two Z bosons to 
decay to electrons or muons is 4.52 x 10~ 3 . The scale 
factor to take into account differences in triggering, re- 
construction and identification efficiencies between data 
and simulation is 0.80±0.08, and the integrated luminos- 
ity is 5.91 ± 0.35 fb _1 . Experimentally, we observe pp — ► 
Z/j* Z/-f* — > £ + £~£ + £~ , and to compare our measure- 
ment with the theoretical prediction of pp — » ZZ, calcu- 
lated in a narrow pole approximation [25], we account 
for Z/^j* interference. The interference in the region 
76 < Mu < 106 GeV/ c 2 increases the acceptance by a 
factor of 1.03. From simulation, the fraction of ZZ events 
that falls outside the region 76 < Mu < 106 GeV/ c 2 is 
0.07 and is also corrected for. The eight observed events 
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therefore result in a cross section: 

a{pp^ZZ) = 2.3 t° Q 9 8 (stat.) ± 0.2(syst.)pb 



where the statistical uncertainty is the 68% confidence 
interval given by the method of Feldman and Cousins 
|26| . The value is consistent with the theoretical pre- 
diction 1.4 ± 0.1 pb [35] • A more precise measurement 
of the ZZ cross section, which combines four-lepton and 
leptons plus fi T channels, is reported elsewhere [27] . 

Examining the properties of the eight ZZ candidate 
events we find an excess of events over standard model 
expectations at high invariant mass, Mzz- The in- 
variant masses of four events are clustered with mean 
327GeV/c 2 , as shown in Fig.JSJ All four candidates, 
one eeee, one ee/z/x, and two n^nn, have values of Mzz 
within 7GeV/c 2 of the mean. In the four-lepton chan- 
nel the detector resolution in Mzz-, c(Mzz)> is 5 to 
6 GeV/ c 2 , so within detector resolution the masses of all 
four events are consistent with a potential new resonance. 
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FIG. 3: M zz for eight ZZ -> l+Tl+T candidates (pythia 
normalized to the standard model prediction of 5.5 events). 



To study the possibility that these events are due to 
a decay of a heavy resonance, we split the eight candi- 
date events into low- and high-mass samples and com- 
pare the properties of the events in the two samples. 
The high-mass region is defined by an a posteriori choice 
Mzz > 300GeV/c 2 , which is ~ 5a(M zz ) below the ob- 
served clustering of events; less than 25% of the expected 
standard model Mzz distribution lies above this cutoff. 

The masses of the Z boson candidates for all events are 
shown in Fig. [4j which demonstrates that the resolution 
in Mu is consistent in the high-mass and low-mass events. 
Lepton identification variables are consistent with expec- 
tation for all the observed events. Most kinematic distri- 
butions for the ZZ — > £ + £~£ + £~ candidates are in agree- 
ment with standard model expectations; as one example, 
the pt distributions of the 16 Z boson candidates are 
shown in Fig. [5] 

However, for the high-mass events, the px distribu- 
tion of the four-lepton system is rather different from the 
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FIG. 4: Invariant masses of dilepton pairs for eight ZZ can- 
didate events: (a) Mu(l) versus Mu(2), with selected mass 
region outlined; and (b) Mu for all Z boson candidates. 
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FIG. 5: pt(Z) for Z boson candidates in (a) low-mass four- 
lepton candidate events and (b) high-mass events (pythia 
prediction normalized to four events in each plot). 



standard model expectation, as shown in Fig. [6] The ZZ 
system in the high-mass events is seen to be boosted and, 
as shown in Fig. [7J is recoiling against one or more jets. 
None of the four low-mass events has a reconstructed jet 
with E T above 20 GeV. 



CD 
o3 



1 r 



CDF, L=6 ft ' 



Mz Z < 300 GeV/c 2 
A data 
^ZZ » 



o4 
> 

CD 
03 



CDF, L=6 It) 1 > 300 GeV/c 2 

▼ data 



(b) 



50 



100 150 200 
p T (ZZ) (GeV/c) 



50 



100 150 200 
p T (ZZ) (GeV/c) 



FIG. 6: pr(ZZ) for (a) low-mass four-lepton candidate events 
and (b) high-mass events (pythia prediction normalized to 
four events in each plot). 



We check whether there is any indication of misrecon- 
struction in these events. In ZZ — > £ + £~£ + £~ events, 
where there is no real fi T , large measured fi T could 
indicate misreconstruction. However the presence of 
jets broadens the detector fi T resolution and needs to 
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be taken into account. To this end we exploit two 
physics models. The first model is RS graviton produc- 
tion through gluon-gluon fusion (the 's-channel signal 
model') [7J. In order to investigate effects of the pro- 
duction mechanism and in the absence of a particular 
model that would predict the production of a boosted 
ZZ resonance, we take as an alternative signal model 
the production of a Kaluza-Klein excitation of a gravi- 
ton, G* , of Mq* — 325GeV/c 2 recoiling against a par- 
ton of Et > 100 GeV (referred to as the 'boosted signal 
model'). In both cases the herwig event generator is 
used with the full CDF detector simulation. In the four- 
lepton decay channel, neither of these models generates 
real fi T . Fig. [7](b) thus demonstrates that the observed 
fi T in the high-mass events is consistent with resolution 
effects arising from the jets. 

Overall, we conclude that the observed events are well- 
measured and that, within the detector resolution, the 
kinematic parameters of the Z candidates are recon- 
structed correctly. The event properties are given in Ta- 
bleU 
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FIG. 7: (a) Number of jets and (inset) Et of the most en- 
ergetic jet; and (b) E T for four-lepton candidate events. E T 
distribution for G*(+jet) process is normalized to 4 events. 



To quantify consistency between the data and the stan- 
dard model, we compute the probability for the observed 
Muu distribution to be due to a statistical fluctuation 
of the standard model expectation. Eight-event pseudo- 
experiments are drawn from the standard model Mzz 
distribution, and a test statistic is computed for each 
pseudoexperiment. 

Two tests are performed. First, the Kolmogorov- 
Smirnov (KS) distance is taken as the test statistic, with 
the intention of testing for goodness-of-fit in a general 
way. The fraction of pseudoexperiments that has KS 
distance greater than that of the observed data distribu- 
tion determines the computed p- value, which is found to 
be 0.14. 

Second, a more powerful test statistic for a resonance 
search is used: the ratio of likelihoods of two hypothe- 
ses. The background hypothesis is provided by the stan- 
dard model distribution in Mzz, Mjfjf, and the signal 
hypothesis adds to it a resonance represented by a Gaus- 
sian peak: / ■ Mff + (1 - /) • G(M,w). For a given 
mass M, the resonance width w is defined by the detec- 



tor resolution at this mass. The resonance parameters are 
defined from fitting the pseudoexperiment distribution in 
Mzz- The likelihood ratio for the data is computed using 
the same procedure. The fraction of pseudoexperiments 
that has likelihood ratio Lsm/Lsm+g lower than that of 
the observed data distribution determines the computed 
p-value and is found to be (1 — 2) x 10 -3 , where the 
range comes from shape differences of the pythia and 
MC@NLO+HERWlG event generators. 

In the absence of a physics model that would predict 
the observed px(ZZ) distribution, we quantify consis- 
tency between the data and the standard model by com- 
puting the probability for eight events sampled from the 
standard model px{ZZ) distribution to have KS distance 
greater than that observed in the data. The probability 
for the data to represent the standard model distribution 
is (1 - 2) x 10~ 4 . 



V. ZZ -» l+Tvv CHANNEL 

The four-lepton events observed above 300 GeV/ c 2 ap- 
pear somewhat anomalous. If these events were to be 
due to a new ZZ resonance, it would also be detectable 
in the other ZZ decay modes, llvv and ££jj. Z bosons 
coming from the decay of such a heavy particle would 
be boosted, so events with one of the Z bosons decaying 
into neutrinos would have large fi T . For each lepton fla- 
vor, the branching ratio to neutrinos is about twice that 
of charged leptons. With all three neutrino flavors in- 
cluded, and only one Z boson to be reconstructed, the ex- 
pected event yield is around ten times higher than in the 
four-lepton channel, and the sensitivity to new physics 
at Mzz = 325 GeV/ c 2 is several times better than in the 
four-lepton channel. 

Optimising sensitivity for a resonance of mass Mzz ~ 
325 GeV/ c 2 , we define the search region to be fi T > 
100 GeV. The standard model expectation for events 
with a Z £ + £~ candidate and such high ft T is around 
25 events, as given in Table |j"l{ Z — s- e + e~ and Z —> 
candidates are selected according to the requirements de- 
scribed for the ZZ — > £ + £~£ + £~ channel. Owing to the 
extra acceptance, we did not reprocess the ££ + JE T data. 

We validate the background model using events with a 
reconstructed Z boson and E T < 100 GeV. Irreducible 
background contributions to a search for new physics in 
this channel come from standard model diboson produc- 
tion processes WW, WZ, and ZZ, as well as from ft 
production. Other non-negligible background contribu- 
tions come from Z+jets events that have large E T due 
to jet mismeasurement; from W+jets events where one 
of the jets is misreconstructed as a lepton and forms a 
Z boson candidate with the charged lepton from the de- 
cay of the W boson; and, in the ee + fi T channel, from 
W'y production with the photon misreconstructed as an 
electron. 

Irreducible backgrounds are estimated using the 
pythia generator and the full CDF detector simulation, 
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TABLE I: Properties of the four-lepton candidate events, in the order in which they were recorded. 



lcptons 


Mz x ,pr{Zi) 


Mz 2 ,Pt(Z 2 ) 


Mzz 


p T (ZZ) 


P T 


Nj e ts 


Jet E T 




(GeV/c 2 ), (GeV/c) 


(GeV/c 2 ), (GeV/c) 


(GeV/c 2 ) 


(GeV/c) 


(GeV) 




(GeV) 


eeee 


93.3, 18.2 


92.9, 17.4 


196.6 


35 


14 







WW 


85.9, 101.9 


92.1, 54.8 


321.1 


47.4 


8.4 


1 


36.7 


ee^i/i 


92.0, 156.0 


89.9, 139.7 


324.7 


126.8 


31 


2 


97.4, 40.0 


eeee 


101.3, 57.8 


91.6, 13.2 


334.4 


44.7 


9.9 


1 


22.7 


eefifi 


87.9, 17.7 


91.8, 29.8 


191.8 


31 


10.5 









95.9, 197.9 


92.0, 87.2 


329.0 


110.9 


23.3 


2 


97.2, 24.7 


eefifi 


95.2, 36.7 


89.7, 38.8 


237.5 


10.2 


1.2 







WW 


88.4, 51.0 


89.8, 26.6 


194.1 


25.9 


3.3 








normalized to NLO cross sections [55]. The Z+jets con- 
tribution is also estimated using pythia simulation and 
is normalized using a subset of the fi T < 100 GeV data. 
As Z+jets events have high fJ T only through misrecon- 
struction, the normalization is carried out on events hav- 
ing 50 < fi T < 100 GeV that also have a small angle 
A</> m in between the fi T and the closest jet, or lepton, 
reconstructed in the event: |A0 min | < 0.5. The |A(/> min | 
distribution is shown in Fig. |8ja). It is verified that this 
procedure is not sensitive to the fi T range used. 
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FIG. 8: (a) A</> m i n as used for Z+jets normalization, and 
(b) Mm for same-sign dielectron pairs with large ]E T used to 
validate the W+jets background estimation. 



ing the same misidentification rates to W ± — > e ± v sim- 
ulation normalized to the NLO production cross section. 
The two methods give results consistent within 10%. 

The overall modeling of the sample composition is 
demonstrated by the fi T spectrum shown in Fig. [9] The 
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FIG. 9: fi T distribution for events with opposite sign lepton 
pairs (ee + /i^t). The contribution of Z+jets events is normal- 
ized in the region 50 < fi T < 100 GeV using events with low 

|A0 m i n |. 



The background contribution from the IF+jets process 
is estimated from a data sample where events contain an 
identified lepton and an additional jet. These events are 
weighted by jet-to-lepton misidentification rates as de- 
scribed in Section IV to estimate the total yield. Owing 



to differences in jet-to-lepton fake rates between electrons 
and muons, the IF+jets contribution is found to be neg- 
ligible in the fifi + fi T channel, but non-negligible in the 
ee + fi T channel. 

Photon conversions are the primary source of jets be- 
ing misidentified as electrons, and so IF+jets events re- 
sult in approximately equal numbers of same-charged 
and oppositely-charged candidate events. The estimate is 
therefore validated against the sample of events that have 
two lepton candidates of the same charge and 50 < fi T < 
100 GeV. Fig.|8jb) shows that this selection is dominated 
by W+jets. The estimate is also cross-checked by apply- 



largest relative uncertainty in this channel comes from 
the Z+jets normalization, and is 10% and 13% in the 
electron and muon channels respectively. Other uncer- 
tainties come from lepton identification (2%), acceptance 
(<1%), cross sections of diboson and top-quark produc- 
tion (5% and 10%), and the fake lepton background 
(20%). The total background uncertainty is 13%. 

To search for a high-mass resonance we examine events 
with fjrp > 100 GeV. Event yields are given in Table |H| 
In electron and muon channels combined we expect 26 
events from standard model processes, and observe 27. 
Four four-lepton events around Mzz = 325 GeV/c 2 com- 
ing from the decay of a new state would imply a produc- 
tion cross section times branching ratio to ZZ close to 
1 pb, and for that cross section, the s-channel G* signal 
model predicts around 35 additional events. 
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As the second Z boson in this channel decays into 
neutrinos, the invariant mass of the Z pair cannot be 
fully reconstructed. The closest approximation is the 
'visible mass' M^f, defined as the invariant mass of 
the sum of the two charged lepton four-momenta and 
the four-vector representing the fi T , (ft x , ft y ,0,\ft T \). 
Fig. 10 shows the M z l z distribution in the signal region, 
.00 GeV, with the expected distribution for an RS 
325 GeV/ c 2 and cross section 
In this channel 



Pt > 

graviton of mass Mq* 
times branching ratio of 1 pb overlaid, 
we find little difference in expected distributions or yields 
between the two signal models, confirming that the anal- 
ysis is not strongly dependent on the detail of the models. 
Neither the event counts of Table [TI] nor the distributions 
of Figure |10| show any evidence for a resonance decaying 
into ZZ. 




FIG. 11: Number of jets in (a) Z — > e + e + > 2jets and 
(b) Z — > fi + fi~ > 2jets events in the control region Mun < 
300 GeV/ c 2 . 
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FIG. 10: Mil for (a) the electron and (b) muon 
channels. The expected contribution from a graviton of 
M G . = 325 GeV/ c 2 and cross section times branching ratio 
to ZZ of 1 pb is shown together with the expected contribu- 
tion of boosted G* , produced in association with a jet. The 
high values of A/Jz of three events in the electron channel are 
understood as originating from fluctuations of the jet energy 
losses in events with high jet activity. 



TABLE 

11+ fi T 



II: Expected and observed event yields in the 
channel. 



Source 


electron channel 


muon channel 


ZZ 


1.8 


1.3 


WZ 


3.G 


2.8 


WW 


0.9 


0.5 


tt 


3.2 


2.4 


W+jets 


0.1 


0.3 


Z+jets 


4.0 


5.1 


Total standard model 


13.6+1.8 


12.4 + 1.6 


Data 


18 


9 


Expected s-channel signal, 
M G = 325 GeV/ c 2 and <r=lpb 


17+1 


18 + 1 


Expected boosted signal, 
Mq = 325 GeV/c 2 and <r=l pb 


20+1 


17+1 



VI. ZZ -> £+r J J CHANNEL 

The decay of a heavy particle into two Z bosons where 
one of the Z bosons decays into charged leptons and 
the other to jets has the advantage of being fully re- 
constructible, and the event yield in the ££jj channel is 
expected to be around twenty times higher than in the 
four-lepton channel. 

Z — > e + e~ and Z — > fi + /i~ candidates are selected 
according to the requirements described for the 
ZZ — > £ + £~£ + l~ channel, and a further requirement is 
made of at least two reconstructed jets having corrected 
i?T>25GeV. To reconstruct the second Z boson can- 
didate, all pairs of jets are considered and if there is a 
pair with invariant mass between 70 and 110 GeV/c 2 it is 
accepted. This inclusive selection, with the additional re- 
quirement of the invariant mass of the two Z candidates 
being less than 300 GeV/ c 2 , defines a control region. 

This channel is dominated by Z+jets events. Other 
standard model sources, small compared with Z+jets, 
are WZ and ZZ production, and tt production. The 
contributions from WW and IF+jets events are negligi- 
ble. 

Diboson and tt event yields are estimated using 
pythia Monte Carlo normalized to NLO cross sections. 
Z+jets events are modeled using the generator ALPGEN 
[2"5] interfaced with pythia for parton showering and 
hadronization, and the normalization of the Z+jets con- 
tribution is obtained by fitting to the total data yield in 
the control region. The detector acceptance is different 
for Z — > e + e~ and Z — > and so the Z+jets normal- 

ization factors for the two channels are not expected to 
be identical. The difference between them is indicative of 
the systematic uncertainty, leading to a total background 
uncertainty of 10%. The jet multiplicity distributions in 
the control region, shown in Fig. |11[ demonstrate the 
good background modeling. 

In the ££jj final state we improve the resolution in the 
reconstructed Mzz by varying jet four-momenta within 
their uncertainties and constraining the reconstructed in- 
variant masses Mjj to the mass of the Z boson, Mz- 
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The resolution in Mz for Z —> jj is 15GeV/c 2 , which 
is much larger than the intrinsic width of the Z boson. 
In the ££jj channel the constraining procedure therefore 
improves the mass resolution of the ZZ candidates, to 
12GeV/c 2 for M G * = 325GeV/c 2 . As the detector reso- 
lution for Mz in Z — > £ + £~ is comparable with the intrin- 
sic width of the Z boson, applying the mass-constraining 
procedure to the leptons has very little effect on the M Z z 
resolution and is used only as a cross-check. Through- 
out this paper Mujj refers to the constrained four-object 
invariant mass. 

To search for a high-mass resonance we examine the 
complete M^jj spectrum. Z bosons coming from the 
decay of a heavy particle would be boosted, and opti- 
mization studies result in requiring the most energetic 
jet in the Z — > jj candidate to have Et > 50 GeV and 
the px of either the Z — > jj or Z — > £ + £~ candidate 
to be greater than 75 GeV/ c. Observed event yields 
are given in Tabic |III| and are consistent with standard 
model expectations. A resonance of Mg. = 325 GeV/ c 2 
and cross section times branching ratio to ZZ of 1 pb 
would be expected to yield around 30 events in the muon 
channel and 40 events in the electron channel, and as 
the ZZ — > £ + £~ jj final state is fully reconstructed, they 
would appear as a narrow peak in Mujj- Fig. 12 shows 



the Mujj distribution for the eejj and fJ.fijj channels, 
with the standard model and additional ZZ resonance 
model predictions. 



TABLE III: 
channel. 



Expected and observed event yields in the Ujj 



Source 


electron channel 


muon channel 


ZZ 


6 


5 


wz 


17 


12 


tt 


7 


5 


Z+jets 


395 


244 


Total standard model 


424±40 


266±24 


Data 


392 


253 


Expected signal, 






M G = 325 GeV/c 2 and <r=l pb 


41±1 


32±1 



Studies of systematic effects resulting from the gener- 
ator Q 2 scale choice and from the jet energy scale uncer- 
tainty show that they do not affect the expected shapes 
of the Mujj distributions. We investigate potential ef- 
fects of the production mechanism using the alternative 
boosted G* signal model. Motivated by the anomalous 
Pt(ZZ) distribution shown by the events in the four- 
lepton channel, the signal selection is modified to require 
Pr(££jj) > 40 GeV/ c, which further suppresses standard 
model background, The resulting Mujj distribution and 
boosted G* prediction is shown in Fig. [l3j As with the 
££ + fi T channel there are no statistically significant dif- 
ferences from the standard model expectation. 
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FIG. 12: Mujj for the (a) electron and (b) muon chan- 
nels, showing the expected contribution from a graviton of 
M G * = 325 GeV/ c 2 and cross section times branching ratio 
to ZZ of 1 pb. 



VII. LIMITS 

To quantify results of the search we compute expected 
and observed limits on the production cross section times 
branching ratio cr(pp —> G* — >• ZZ). 

The expected sensitivity is determined with a Bayesian 
technique using CLg likelihood test statistics [3D] 
to perform a binned maximum-likelihood fit over the 
Mzz, Mzzi an d Mujj distributions in the four-lepton, 
££ + fi T , and Ujj channels respectively. The background 
hypothesis is provided by the standard model expectation 
as described in Sections |IV|VI Background-only pseudo- 
experiments are drawn from Monte Carlo simulation. In 
the fit, the background templates can fluctuate within 
their uncertainties. A test statistic is formed from the 
difference in the likelihoods between the background-only 
model and the signal-plus-background model at the best 
fit values for the pseudoexperiment. From this, expected 
95% credibility level (CL) upper limits on cross section 
times branching ratio are extracted. 

Fig. [14] shows expected and observed limits in the 
four-lepton channel for G* masses between 250 and 
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FIG. 13: Mujj for the (a) electron and (b) muon channels for 
Pt(ZZ) > 40GeV/c, showing the expected contribution from 
a boosted graviton of Mg* = 325 GeV/ c 2 and cross section 
times branching ratio to ZZ of 1 pb. 



Mg* = 325GeV/c 2 , and the observed limit is 0.23 pb. 
With the selection tuned for a boosted signal model, 
Pri^jj) > 40GeV/c, the sensitivity is improved slightly 
compared to the s-channel signal model. The expected 
limit is 0.27 pb and the observed limit is 0.26 pb, show- 
ing that also in this channel the analysis sensitivity is not 
strongly dependent on the detail of the signal model. 
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1000GeV/c 2 . At M G * = 325 GeV/ c 2 the expected sen- 
sitivity is around 0.7 pb, and the four events with masses 
clustered around that value result in an observed limit of 
1.9 pb. 



41 expected i 

□ 1-a 

□ 2-a 

— 41 observed : 
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FIG. 14: Expected and observed 95% CL limits on <j(pp — > 
G* — > ZZ) from the four-lepton channel; the four events with 
Mzz=327 GeV/c 2 result in a deviation from the expected 
limit. 



Although the backgrounds in the 11 + ft T channel are 
higher than in the four-lepton channel, this channel pro- 
vides better sensitivity. Fig. 15 'a) shows the expected 
and observed cross section limits for 11 + fi T , and there 
are no large differences from standard model expecta- 
tions. For Mg* = 325GeV/c 2 the expected 95% CL up- 
per cross section limit on the s-channel signal model 
is 0.29 pb, and the observed limit is 0.25 pb. For the 
boosted G* signal model the 95% CL expected and ob- 
served limits are both 0.30 pb. This is a change of less 
than 10% from the s-channel model, demonstrating that 
the analysis sensitivity is not strongly dependent on the 
detail of the production model. 

Fig. 151)) shows the expected and 
limits for the Iljj channel. 



observed cross 
Here the ex- 



FIG. 15: Expected and observed 95% CL limits on er(pp — > 
G* -> ZZ) from (a) the ZZ -¥ l + l~vv channel, and (b) the 
ZZ -> £ + r jj channel. 



Combining all three channels results in the most sensi- 
tivity. Expected and observ ed l imits are consistent with 
each other, as shown in Fig. 16 For Mg- — 325 GeV/ c 2 



the sensitivity is dominated by the 11 + fi T channel. For 
an s-channel resonance, the 95% CL upper cross section 
limit is expected to be 0.19 pb and is observed to be 
0.26 pb. For a boosted resonance of Mg. = 325GeV/c 2 
the expected limit is 0.17 pb and the observed limit is 
0.28 pb. 
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FIG. 16: Expected and observed 95% CL limits on cr(pp 
G* — > ZZ) from all channels combined. 



section 

pected 95% CL upper cross section limit is 0.38 pb for 



VIII. CONCLUSIONS 

We have searched for heavy resonances decaying into 
Z boson pairs using the final states consisting of four 
leptons, two leptons and ft T , and two leptons plus jets. 
In the channel with the smallest background, the four- 
lepton channel, we have observed eight candidate events. 
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Four events with high values of ZZ mass are close in 
mass, and two of those have unusually high pt(ZZ). 

However, more sensitive searches in the it + fi T and 
ttjj final states show no indication of a new heavy parti- 
cle decaying to two Z bosons, suggesting that the events 
observed around 325 GeV/ c 2 in the four-lepton chan- 
nel result from standard model processes. Combining 
all three channels we set upper limits on the cross sec- 
tion times branching ratio a(pp — > G* — > ZZ) that 
vary between 0.26 pb and 0.045 pb in the mass range 
300 < M G , < 1000GcV/c 2 , and the limits do not de- 
pend strongly on the production model. 
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